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P_H | Abstract 

' We calculate gluon production in deep inelastic scattering of the current j = —^F^ U F^ V off a large 

^Nj . nucleus and in nucleon-nucleus collisions. In a covariant gauge calculation the transverse momentum 

spectrum of the gluon is determined by the final state interactions of the gluon with the nucleons in the 

nucleus. In a light-cone gauge calculation final state interactions are absent and these effects come from 

the light-cone wavefunction of the nucleus. We work in an approximation which neglects QCD evolution 

^-T ' of gluons in the nucleon, that is in a quasi-classical approximation. 

'xT ■ 

-xf ; 

^ ■ 1 Introduction 

00 ■ __ 

A few years ago L. McLerran and R. VenugopalanpJ introduced an interesting program of studying the 

small x gluon distribution, and gluon production, involving a very large nucleus. The essential idea is to 
use the valence quarks of the nucleus as a source for a light-cone gauge calculation, involving only tree 
graphs, of the Weizsacker- Williams (WW) field of the nucleus. The higher momentum components of this 
WW field, interpreted in terms of gluon quanta, would then be included with the valence quarks to give a 
new source which could be used to calculate another momentum layer of WW quanta, etc. This program 
is being actively pursued [|| [|. The calculation of the WW field of the valence quarks of a large nucleus is 
. | a nontrivial problem because of the nonlinearity due to the strong source coming from the many valence 

■ quarks of the large nucleus. It was calculated independently in || and ||. 

It was observed in this (quasi-classical) WW approximation that the gluon distribution for a large nucleus 
is strongly modified from what one would expect from adding the distributions of the individual nucleons 
incoherently. An apparent saturation of the gluon distribution sets in at moderate to low gluon transverse 
momenta while the high transverse momentum part of the WW single particle gluon distribution remains 
additive. However, some time ago|| it was observed that there are no shadowing corrections, in this 
equivalent gluon approximation. It is one of the purposes of this paper to show the compatibility of 
the no shadowing result with the result of Ref.5 and to interpret that latter result. To do this we introduce 
a current j(x) = — \F^ v (x)F^ l/ (x) and calculate deep inelastic scattering of this current off a large nucleus 
in an approximation where the nucleons in the nucleus have no QCD evolution included in their gluon 
distributions. We suppose a <C 1 but the parameter a 2 A 1 / 3 , with A the atomic number, will generally be 
taken to be large. We do our calculation first in covariant gauge and then in an appropriate light-cone gauge 
in order to interpret the result of Ref.5. 

It is convenient to consider two distinct regimes of gluon production in the deep inelastic scattering of 
the current j off a large nucleus. If Q 2 = —q^Qfi is the virtuality of the current and if a/< > is the typical 
transverse momentum that a high energy gluon obtains by multiple scattering with the nucleons as it passes 
through the nucleus H, then these two regimes are characterized by Q 2 ^$>< l\ > and by Q 2 < < £ \ > 

2 
J_ 
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We begin by describing our results when Q 2 ^>< i\ > 



*This work is supported in part by the Department of Energy under GRANT DE-FG02-94ER-40819. 



1 



In small x deep inelastic scattering, and in covariant gauge, the current j produces a gluon locally in 
the nucle us p| . This gluon then multiply scatters, both elastically and inelastically, as it passes through the 
nucleus^, |9|, |l0|. The equation for the transverse momentum distribution of the produced gluon is given in (6) 
below. The solution, given in (18) and (21), is exactly the one found in || for the gluon correlation function. 
This distribution does not reflect shadowing, but rather the probability conserving final state interactions 
which modify the transverse momentum distribution of the produced gluon but not the cross section for its 
production. In an A + = light-cone gauge calculation, in a frame where the nucleus is right-moving and the 
current left-moving, we find a complete absence of final state interactions and we give a physical picture of 
why this happens. Here the effect of the final state interaction which we found in covariant gauge is encoded 
in the light-cone wavefunction of the nucleus. It is quite remarkable that the distribution of gluons in the 
WW wavefunction gives the correct transverse momentum distribution of produced gluons without any final 
state interactions whatever. The traditional gluon distribution xGa(x, Q 2 ) can be obtained by integrating 
the produced gluon distribution over all transverse momenta obeying i 2 ^ < Q 2 . Because Q 2 3> in the 
logarithmic approximation, one can view the current j as measuring the gluon distribution in the nuclear 
wavefunction. 

When Q 2 < < t\ > the covariant gauge calculation of the produced gluon spectrum does not change in 
form from that when Q 2 ^$>< l\ > since the large < £\ > comes from final state interactions which appear 
after the production of the gluon occurs. The light-cone gauge calculation now appears more complicated. 
The three and four gluon terms in j become important and final state interactions are not negligible. 
The current is no longer sufficiently pointlike to be viewed as measuring the partonic gluon distribution of 
the nucleus. If we construct xGa(x,Q 2 ) from its moments, defined in terms of nuclear matrix elements 
of the normal local gluon operators, with the counter terms taken to be the same as for nucleon matrix 
elements, then xGa(x, Q 2 ) is additive xGa(x,Q 2 ) = AxG(x,Q 2 ). It appears that this object has no direct 
connection with partonic gluons in the light-cone wavefunction of the nucleus when Q 2 < < i\ > . It is 
perhaps surprising that the gluon distribution in the light-cone wavefunction is still the same as the produced 
gluon spectrum. Finally, while the gluon distribution in the light-cone wavefunction is very different from 
the incoherent nucleon light-cone distributions when i\ < < i\ > the nuclear modifications are simply 
probability conserving redistributions in phase space and so are unrelated to nuclear shadowing. 

In Sec. 3 we calculate gluon production in a nucleon-nucleus collision. In a calculation equivalent to a 
covariant gauge calculation we view the gluon as being produced in one of two ways, (i) The gluon may 
be present in the wavefunction of the hadron as the hadron reaches the front face of the nucleus. In that 
case interactions with nucleons in the nucleus free the gluon and multiple scattering in the nucleus broadens 
its transverse momentum spectrum, (ii) The gluon may be radiated in the final state from fragments of 
the incident hadron. In that case the transverse momentum spectrum of the gluon is not affected by the 
nucleus. Interference terms between initial and final state emission are small. In an A + — light-cone gauge 
calculation of gluon production from the collision of a right-moving nucleus with a left-moving nucleon no 
final state interactions appear. If is the field of the nucleus while A 1 ^ is the field of the nucleon gauge 
rotatedj 11 by those nucleons in the nucleus having the same impact parameter then the gluon production 
cross section is obtained by coupling A^ and A'^ to the QCD radiation field in a way determined by the 
QCD Lagrangian. A very simple expression emerges for the cross section as given by (73)- (75). 



2 The gluon distribution and gluon production in deep inelastic 
scattering off a large nucleus 

In this section, we consider deep inelastic scattering off a large nucleus where the (gauge invariant) current 
initiating the scattering is taken to be j(x) = — \F^ V F^ V with F^ v the usual QCD field strength tensor. 
We choose this current because it furnishes the most direct way to measure the gluon distribution and to 
produce a gluon and follow its interactions with the nucleons of the nucleus as it passes into the final state. 
We neglect QCD evolution in the interaction of j with the individual nucleons in the nucleus. Our calculation 
of the gluon distribution is closely related to that of Ref.5 with which we are in agreement. Our motivation 
for neglecting QCD evolution, and we shall explain below in what sense we neglect it, is the same as in Ref.l. 
We are interested in a model to illustrate particular small x and nuclear effects. We shall later comment on 
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our expectations in a realistic nucleus. 

There is a surprisingly different picture of gluon production in covariant gauge as compared to light-cone 
gauge. In covariant gauge, and in the rest systems of the nucleus, the current j produces a gluon off a nucleon 
in the nucleus in a reasonably local way. This gluon then has final state interactions with nucleons along its 
path as it passes through the nucleus jlO|, |l2). The final state interactions change the transverse momentum 
distribution of the produced gluon, but they do not modify the production rate. Thus the gluon distribution, 
and the structure function, of the nucleus is additive in the atomic number and there is no shadowing. In 
light-cone gauge, at least in a particular light-cone gauge, there are no final state interactions whatever. 
The transverse momentum distribution of the produced gluon is already encoded in the wavefunction of the 
nucleus. After deriving this result we shall give a simple physical picture why final state interactions are not 
present in our light-cone gauge. 

2.1 Covariant gauge calculation of gluon production of the nuclear gluon dis- 
tribution 

We begin with the scattering of j = — \F 2 V off a single nucleon as illustrated in Fig.l. Define a structure 
function 

W = (2n) 3 2E p J d 4 x e**(p\j(x)j(0)\p) (1) 

in analogy with deep inelastic lepton nucleon scattering. Then, from the graph in Fig.l it is straightforward 
to show 



W(s,Q 2 ) = -W £n s/p 2 (2a) 



which can be factorizcd as 



W{s,Q 2 ) = -xG(x,Q 2 ) + -W £n 1/x (2b) 

with x the usual Bjorken variable, s = (p + q) 2 and Wq a constant. (The Feynman rules for coupling j to 
two gluons are given in Ref.6.) The first term on the right-hand side of (2b) corresponds to gluon operators 
in the operator product expansion while the second term corresponds to quark operators. As discussed in 
some detail in Ref.6 the nucleons are additive in our approximation so that 

W A = AW (3) 

and 

xG A = AxG (4) 

Where Wa and Ga are the structure function and the gluon distribution for the nucleus. Indeed, additivity 
is closely related to our treatment of the gluon- hadron amplitude, the lower blob in Fig.l. We suppose 
that there are no £n k 2 /A 2 factors in that amplitude and we also assume that the gluon-nucleon amplitude 
vanishes as k ■ p becomes large. That is, we suppose that high mass intermediate states are not present 
in the gluon-nucleon scattering amplitude. These assumptions limit the longitudinal coherence length of 
the forward current-nucleon scattering amplitude to be not significantly larger than the nucleon size, which 
means that the interaction is local in the nucleus. 

Now let us follow a gluon, produced by an interaction of j with a nucleon in the nucleus, as it scatters 
with other nucleons in the nucleus as illustrated Fig. 2 where we have shown examples of both inelastic and 
elastic scattering in the nucleus. It is convenient to define a probability distribution for the gluon to have 
a transverse momentum I as it passes through the nucleus JT^. If the gluon is produced with a transverse 
momentum £ at a longitudinal coordinate zq and at an impact parameter b with respect to the center of 
the nucleus, let f(b, zo,£ , z,£) be the probability distribution for the gluon to have transverse momentum £ 
at a longitudinal position z. Of course 
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d 2 £f(b,z Q ,£ ,z,£) = l (5) 
while the number density of produced gluons in the nucleus is 

= / d2e o^f/bdz p(b, z )f(b_, z 0l £ , z,£)\ z=VW ^ (6) 

where the integration d 2 bdz goes over the volume of the nucleus, p is the nuclear density normalized to the 
atomic number A by 

d 2 bdzop(b, zq) = A, (7) 

and we take the nuclear density to be uniform throughout the nucleus for simplicity. The initial distribution 
resulting from an interaction of the current j with a nucleon in the nucleus is normalized to be given by the 
unintegrated gluon distribution of the nucleon 

dN (Q ldN Id 2 

-PT ■= = ~w [xG{xA )] (8) 

where, in our approximation, there is no x-dependence in xG(x,£ 2 ) and there is no ^—dependence in 
£ 2 -J^xG{x,£ 2 ) for large i? 2 . 

Now f obeys the equation (See Eqs.2.7 and 2.8 of Ref.12. Appendix A of that reference gives a derivation.) 

j^/M)= -\f{z,Q+pa j d 2 £'V{£!)f{z,£-0 (9) 

with 

f(zo,£) = S(£_-£ ) (10) 

and where 

l '< f > - h% <"> 

is the normalized gluon-nucleon scattering amplitude with £ the momentum transfer and where 

A=H" 1 , (12) 

is the mean free path of gluons in the nuclear matter. For notational simplicity we suppress the b,zo,£ 
dependence of f. The first (inelastic) interaction of the produced gluon shown in Fig. 2 is contained in the 
second term in (9) while the second (elastic) interaction comes from the first term in (9). We work in an 
approximation where the gluon-nucleon cross section is energy independent, at high energy, and where the 
gluon's helicity is conserved during the scatterings in the nucleus. 
Eq.9 is solved by going to a transverse coordinate representation 

f(z,x)= J d 2 £e-W--^f{ z ,£). (13) 

As usual in the high energy limit the transverse coordinate of the gluon will not change as the gluon passes 
through the nucleus, and thus, using (13) in (9) gives 

§ ; f\z,x) = -^x 2 v(x)f(z,x) (14) 

with 

f(zo,x) = 1 (15) 
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and where 



with 



Now defining 



v(x) = ^(l-V{x)) (16) 



V(x) = / d*le-*zV(£. (17) 



N(3L) = / St e-^-^- (18) 

one easily finds from (6) that 

N(x) = f d 2 bdz p(b,z Q )NoU)f(z,x)\ z=vw ^ (19) 



The solution to (14) and (15) is clearly 



(Z-Zo) 2, 



leading to 



f(z,x) = exp{-^^x 2 v(x)} (20) 



N 2 - 1 J R 2 — h 2 

N{ - ] = I d2 K 2 aN x 2 (1 ~ GXP[ 2X (21a) 



,>, N 2 -l ,„ r 27rVi? 2 - b 2 x 2 aN c 



N(x) = / d 2 b 2 C 2 (1 - exp[ ~ ^pxG(x,x 2 )]) (216) 

J ■n^aNcxf N£ — 1 

where we have used (See Eq.4.7 and Appendix B of Ref.12.) 

v(x) 4ir 2 aN c 9 . . . 

valid in a logarithmic approximation for small x 2 . (Our notation is such that xG(x,x 2 ) is not the Fourier 
transform of xG, but rather stands for xG(x, 1/x 2 ).) 

Eq.(21) is the same result found in Ref.5. In the limit of very small x 2 one can keep only the first two 
terms in the expansion of the exponential in (21) to arrive at 

N(x) = AxG(x,x 2 ) (23) 
which, with the identification of the nuclear gluon distribution with N(x\), gives 

xG A (x,x 2 ) = AxG{x,x 2 ). (24) 
At first sight it would seem that (24) no longer holds for extremely large nuclei at fixed x 2 where 

N(x) « tL. (25) 

However (24) is still correct. What fails is the identification of N(x) with xGa(x,x 2 ) when A is very large. 
The physics here is quite clear. The current j makes a (local) hard interaction with one of the nucleons in 
the nucleus and produces a gluon. The coherence length for this production is small and so the production 
cross section, and hence the nuclear gluon distribution xGa(x, Q 2 ), is incoherent among the nucleons in the 
nucleus [p"3"[ . The quantity which is calculated in (6) and (21) is the momentum distribution of produced 
gluons, and the x in (21) is conjugate to the £ of the produced gluon but has no natural connection with 
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the Q 2 of the current when A is very large. When A is not too large, and the final state interactions do not 
change the momentum of the produced gluon too much, and if one takes a cutoff on (? equal to the Q 2 of 
the current, then N(x 2 = A/Q 2 ) can be identified with xG^fi,! 2 ) and (23) holds. However, when the final 
state interactions give a broadening of the transverse momentum distribution A£ 2 which is large compared 
to Q 2 then N(gf = 4/Q 2 ) has no relation to the nuclear gluon distribution. 

2.2 Light-cone gauge calculation of gluon production and of the nuclear gluon 
distribution 

In a light-cone gauge calculation of gluon production one must also arrive at (21) and (24). Indeed (21) 
has already been obtained in Ref.5 using a light-cone gauge although the interpretation of what exactly was 
calculated is, perhaps, not so straightforward. In this section, we calculate gluon production in rj-A = A + = 
light-cone gauge. We restrict our discussion to Q 2 3> £ 2 . We shall find (21) although we shall see that there 
are no final state interactions when I 2 < Q 2 . The transverse momentum distribution of the produced gluon 
now comes from the wavefunction of the nucleus. 

We begin by showing that final state interactions are absent. The only unusual aspect of our calculation 
is the way we choose the it's in the light-cone gauge propagator. With the momentum as indicated in Fig. 3 
our light-cone gauge propagator is 

al3 k 2 + it a/3 rj ■ k — it rj ■ k + it 

This choice was used in Refill and implicitly in Rcfs.4 and 5. If r\ ■ k = k+ then the choice of it's in (26) 
dictates that a charge at an x_ = x "^ 3 coordinate gauge rotates all charges having X- < x^ but does 

not gauge rotate charges with x_ > x^ . 



2.2.1 Final state interactions are absent 

In order to understand why final state interactions do not appear in deep inelastic scattering in light-cone 
gauge it is useful to begin with the example of gluon production off a quark in one nucleon and the subsequent 
rescattering of the gluon off a quark in another nucleon as illustrated in Fig. 4. The quarks are assumed to 
be right-movers and the produced gluon, £, is assumed to obey £ + <C The quarks are on-shell before 
and after the scattering while the quark having p2 is taken to have its ^--coordinate to be greater than 
the ^--coordinate of pi so that it is natural to expect a production of the gluon off quark 1 and a later 
rescattering off quark 2. The + symbols along the quark lines in Fig. 4 indicate that the 7+ matrices will be 
dominant, as usual when considering soft gluon production. 
The relevant factors for the graph shown in Fig. 4a are 

°- = ( t + 9 -OT + t ' + ",-V + «l ' ( * - ' + «>»<' - l) " - «"<* - l) ■ {k + " - t)] 



2t_[(k -£) + + ^- it] ( fc - ')+ + ^ ( k l )+ 

—ik 1 - 

r 7 „ g .. e„(l). (27) 

k (fc+ - it) 

We now carry out an integration over fc + , recalling that £ + <C —q+ <C p+.{t+ <C —q+ follows from £- = q- 
and Q 2 = —2q + q_ 3> £? while — q + <C p+ follows from the fact that we are considering small-x deep 
inelastic scattering with x = —q + /p + . We take transverse components of g M to be zero.) Before doing the 
k + integration we use the Slavnov- Taylor- Ward (STW) identities to drop the rjp(k — £) 1 part of the gluon 
propagator as this term cancels between the three graphs, (a-c), of Fig. 4. After this term is dropped only 
the ( fc+g _^ ++ie has a singularity in the lower half fc + — plane. We can do the dk + integration by distorting 
the contour in the lower half k + — plane so long as the rest of the integrand, the other terms in G a , vanish 



G 



for large values of k + . And if that is the case the net result will be small because the final evaluation is done 
at k + = (£—q) + w — q + so that inverse powers of k + correspond to higher inverse powers of Q 2 as compared 
to the case of no rescattering. In order to count the k + — factors in G a it is convenient to write 



9 * ~ (k - £) + Z £ = + *»> + • (28) 



fk— i) 2 

Then in order to compensate for the k + — ie and (k — £)+ + 2l ie denominators, not including the 

possible denominator in (28), one would need a k\ coming from the numerator terms in order not to have 
a small result. This means that one must get a single factor of fc+ (or q + ) each from the T iav — term and 
from the []/3 Q — term in (27). In order to get a k + — factor from r 7<T „ an accompanying fj-y , fj a or fj v must 
appear. But fj a clearly gives zero from (27) while rfy gives zero, or an extra fc+ in the denominator, from 
(28) so we are left with fj v as the only possibility. Here we can simplify our discussion by choosing fj ■ e = 0. 
This is perhaps an unusual choice to go with n ■ A = light-cone gauge but, after all, one is free to choose 
the external polarizations in any convenient basis. With the choice fj ■ e — we find 



/ 



dk+Ga = (29) 



(We could arrive at this result somewhat more directly. After observing that there is sufficient convergence 
to do the k + — integration by distorting in the lower half plane we could use current conservation in the 
a- index of the q— vertex along with (k + q — £)■£ w (k + q — £) a to get (29). However, it does not seem 
possible to generalize this procedure to all higher order graphs.) 

The choice fj ■ e = also makes graphs b and c of Fig. 4 equal to zero because one cannot have a 7+ 
matrix with this polarization vector. Thus we are left only with the graph shown in Fig.4d and we must 
also evaluate it with fj ■ e = 0. Now 

r = -i{k-£ + g)j -ik± 

d {l-l + <\) 2 [{k + q-£) + + ie]\, 2 {k + -ie) lJaav ( ' 

where 



r-U = 9co{£ -q- 2k) u - q au {2£ - q -k)„ + g va {k + £) a . (31) 

Again, in doing the k + integration we will end up with a q + in the denominator, giving a small result, unless 
one Cctii tcike ci k-\~. or a q + , from T as given in (31). But, taking a k + , or a q + , from T will result in a factor 
of fj v , fja or fj a each of which give zero in (30). 

We have seen that the graphs of Fig. 4 give a small contribution. (We have reached the conclusion for 
£ 2 <C Q 2 - Of course we would have found a small result for large independently of the size of Q 2 , but 
that is no different than in covariant gauge. The fact that final state interactions are suppressed for large 
Q 2 , and with if <§C Q 2 , is different than what we found in covariant gauge.) In a moment we shall give a 
physical argument explaining why this is not unexpected. However, before describing our physical argument 
consider the high energy gluon quark scattering amplitude at order a in A + = light-cone gauge and in 
a frame where the quark is right-moving and the gluon left-moving. The graphs are shown in Fig. 5. If we 
choose fj ■ e — for the initial and final gluons, then only the first graph in Fig. 5 contributes. The essential 
factors in this graph are 

G = u(p k) 1+ u(p)^ ^r^e^e^. (32) 

The on-shell conditions {£ — k) 2 = and (p — k) 2 = give k + — ~ 1 = ~ k 2 /2p+- In order to see 

how the high energy limit comes about it is useful to write 

ka = k a - fj a k + - n a k- . (33) 
k a gives zero while acting in (32) while rj a gives a small factor. Using 
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(34) 



gives 



i„ 2 - - 



(35) 



Of course, the result emerges directly from an evaluation of (32) without explicitly invoking current conser- 
vation, (33). In a frame where £- is large the high energy growth, the factor of s in (35), comes partly from 
the smallness of fc+ in the denominator in (32). In order to get a small denominator it is crucial to have 
a long region of integration over X-, the variable conjugate to Of course the quark and gluon are far 
apart when X- is large, but in axial gauge important interactions can happen in a non-causal manner. In 
our choice of light-cone gauge, (26), the factor ^- in (32) comes from times long before the quark and gluon 
have reached each other. 

We are now in a position to see why final state interactions cannot be important in single gluon production 
in light-cone gauge in deep inelastic scattering. In Fig. 4a and its higher order counterparts a final state 
interaction occurs much like the one we have just considered. In order for this final state interaction to be 
important one would need a trapping of the fc+ — integration contour at values of fc+ on the order o{£ 2 /2£_ 

but from (27) we see that we can distort the fc + -contour so that fc+ is always of size q + = Because Q 2 

is a large parameter fc + is not forced to be small, and thus final state interactions are weak. We can say this 
slightly differently. The highly off-shell current has a lifetime Aa;_ ss ^fc- and this is the limiting X- over 
which integrations can be done. This Ax- is not large enough to generate the small light-cone denominators 
necessary to have final state interactions, at least so long as the target is not too long compared to its 
longitudinal momentum, that is so long as £ 2 «Q 2 . 

2.2.2 Gluon production in deep inelastic scattering 

We begin by considering deep inelastic scattering in a nucleus where two separate nucleons are involved in 
the scattering. The two classes of graphs are illustrated in Fig. 6. The graphs of Fig. 6a correspond to inelastic 
scattering off both nucleon 1 and off nucleon 2 while graphs in the class of Fig. 6b correspond to an inelastic 
reaction with nucleon 1 and an elastic reaction with nucleon 2. We emphasize that in our frame where the 
nucleons are right-movers and the current q has q_ 3> q+ the X-- coordinate of nucleon 1 is less than the 
X-— coordinate of nucleon 2. The interactions shown in Fig. 6 are nonzero only when the propagator of the 
fc— line is taken to be _ ie ] and similarly for the (k — A)— line. (We suppose A± = A_ = while 

A + is integrated freely over values large compared to £ + as is appropriate for a low x collision where the 
produced gluon has a coherence length large compared to the length of the nucleus.) Our gauge choice, as 
indicated in (26), allows the gluon field of a nucleon to have a pure gauge part which extends in the negative 
X- direction starting from that nucleon. As shown in detail in Ref.ll this field gauge rotates the field of 
nucleons having smaller values of X-. Thus, the field of nucleon 2 rotates the field coming from nucleon 1. 

Now let us evaluate explicitly the graphs in Fig. 6 starting with those in Fig. 6a. The fc— line of that figure 
goes from nucleon 2 to nucleon 1 or to the field coming from nucleon 1. Because of the k a at the end of the 
fc— line we can use the STW identities to evaluate these contributions. This has already been done in Ref.ll 
and here we need only outline the procedure and put in the normalizations appropriate for our problem. In 
doing the fc+ — integration we distort over the k ^_ ie pole if the light-cone propagator which sets fc+ = 0. 
Since fc_ = compared to £_ the fc— line, which by the STW identities eliminates the [l — g)-propagator, 
effectively inserts a transverse momentum fc just before the attachment of the current q. In terms of formulas 



where T aap is the triple gluon vertex, D pfl the gluon propagator and the vertex of the current q with the 
gluon lines £ — q and t. The omitted terms in (36) cancel with gluon, k a , attachments into the lower blob of 
nucleon 1. In exactly the same way, and integrating over A + , the gluon line fc — A also eliminates the t — q 
line in the complex conjugate amplitude. 




(36) 
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Thus, the STW identities allow one to calculate the effect of nucleon 2 on the field coming from nucleon 
1. We may write the gluon production amplitude in terms of the scattering of the current off the gluon 
(£ — q), as illustrated in Fig. 7, times the unintegrated gluon distribution coming from nucleon 1 as modified 
by nucleon 2. Rather than immediately writing the answer for the two nucleon case we instead proceed to 
the general case using the "classical" field calculated in Rcf.4 to generate the tree graphs contributing to 
the unintegrated gluon distribution coming from an arbitrary number of nucleons in the nucleus. We change 

notation slightly from Ref.4 to better match the problem at hand. Write A^(x, X-) = X^i" 1 ^Af^-fjc, x -) 
in light-cone gauge as 



with 



A^{x,x-) = j S{x,b-)T a S-\x,b-)\7^1n[\x-b\fj]p a (b,b-)d(b--x-)d 2 bdb- 



S(x,x_) = Pcxp^igT a J ln[\x -b\fj]p a (b,b-)9(b- - x_)d 2 bdb- 



(37) 



(38) 



where P path orders the integration with terms having smaller values of 6_ coming more to the right. 
p a is a color charge operator normalized according to 

<p a (b,b.)p a ' (b,'b'_) >= ^M S (b-b')6(b.-b'_)6 aa ,Q 2 ^xG(x,Q 2 ) (39) 
with p(b, b ) the normal nuclear density, in our boosted frame, obeying 

d 2 bdb-p(b, b-) = A. (40) 



/i in (37) and (38) is an infrared cutoff which will disappear from physical quantities. We note that 
Q 2 g^2-xG(x, Q 2 ) is independent of x and Q 2 in our approximation. The expectation indicated in (39) 
is an expectation in the nuclear wavefunction and we suppose that the pairwise correlation indicated there 
is the only nontrivial correlation. This is a natural assumption for a nucleus having weak nucleon-nucleon 
correlations. Define a Fourier transformed field by 

A(k,k+) = J ^ J dx-e i{ - k +- ie >— ik ^A{x,x-). (41) 
Then the main formula which we shall need in order to evaluate ^ is 

tt^ = -2Tr < A^(£, l + q + )Aj;{-l, -£ + + q + ) > [g ^ ~ i? i^Q 2 }^ (42) 



where the term in the brackets on the right-hand side of (42) is given by the graph in Fig. 7. The A^ in 
(42) is the same as given in (37) and represents the field of the nucleus which the source j interacts with. 
Eq.(42) does not have any final state interactions, interactions of the produced gluon with the nucleus after 
the action of j, which we have previously argued are small. As Q 2 ^> £ 2 and l + <C — q+ then taking the 
leading Q 2 term in Eq.(42) and using Eq.(18) yields 

N(x 2 ) = d2bTr < A^(b)A^(b + x) > . (43) 

Plugging in Aj^ from Eq.(37) one finds 

N(x 2 ) = ~J d 2 b j d 2 b'db'_d 2 b"db"_ 
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Tr[S{b,b'_)T a S- 1 (b,b'_)S{b + x,b"_)T a S- 1 (b + x,b'L)} > . (44) 

Here we can assume that the path ordered 6_ -integration in the S(x, x_)-matrices [see Eq.(38)] is performed 
up to some point close to x~, but excludes this point and its immediate vicinity. This statement is equivalent 
to what was referred to as dropping the last nucleon in Ref.4. Using these assumptions we can independently 
calculate the correlation of the two color densities in (44). Employing (39) and integrating over cPb' we arrive 
at 

N(x 2 ) = 4 J ff±_ln(\x\ f ,)Q 2 ^ X G(x, Q 2 ) 

J d 2 bdb'_ < Tr[S(b,b'_)T a S- l (b,b'_)S(b + x,b'_)T a S- 1 {b + x,b'_)] >, (45) 

where p, is some infrared cutoff, and wc have assumed that the normal nuclear density is uniform throughout 
the nucleus and in the boosted frame is given by p re i- Rewriting (38) according to the definition of a path- 
ordered product we obtain 

S(b, &_) = + igT a p a {y_, y-)£n(\b - y_\p)d 2 yAy^ 

i 

-(l/2)g*T a T b p a (y, y.)p\y^ y'_)ln{\b ~ y_Wn{\b ~ y » 

d 2 yAy^d 2 y'Ay' l _], (46) 

where the product goes along the y~— axis, such that smaller y_ corresponds to greater label i. For a 
combination of matrices like S" -1 ^, b'_)S(b + x, b'_) we start by considering the last (rightmost) term of the 
product given by (46) for S 1-1 ^, b'_) and the first (leftmost) term in a similar product for S(b + x, b'_). Both 
terms include the same interval of the y_ integration, the closest to b'_. Since all S— matrices in (45) are 
taken at the same longitudinal coordinate b'_, then taking these terms for both pairs of S— matrices in (45) 
we observe that nothing else in (45) depends on the longitudinal coordinates in this interval. Therefore we 
can average these terms independently of the rest of the expression in (45), as well as do the transverse 
integrations. Taking only the two-density correlation terms (We throw away the higher order correlations 
since they have more powers of g, which goes beyond the classical approximation jll]]) and making use of 
(39) we end up with 

< Tr[S(b,b , _)T a S- 1 (b,b , _)S(b + x,b'_)T a S- 1 (b + x,b'_)} > 

< Tr[S{b, b'_ - Ay^S^ib, b'_ - Ay_)S(b + x, b'_ - Ay^S^ik + x, b'_ - Ay-)} >, (47) 

where Ay_ is the absolute value of the change in longitudinal coordinate. The transverse integration d 2 y was 
done assuming that the nucleus is infinite in the transverse direction, which is a reasonable approximation 
for a large nucleus. In arriving at (47) we have neglected QCD evolution in xG. Continuing this procedure 
of picking small intervals along the y_— axis and taking the limit of Ay_ — * shows that the trace in (47) 
is equal to an exponential function 

< Tr[S(b, bf_)T a S- 1 fe b'_)S(b + x, t/_)T a S~ 1 (^ + x, b'_)\ >= 

C F N c exp [-9 2 ^f^f G(x, l/x 2 )(b'_ + b' _f) , (48) 

with ±6q_ the upper (lower) limit of the y_ integration in the S— matrices in (47). Plugging this back in (45), 
performing the db'_ integration and defining the nuclear density in the center of mass frame p = p re i / 7v2> 
with 7 the Lorentz contraction factor provides us with 
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N(2L 2 ) = / d 2 b 

-K z aN r x z 



1 - exp ( - 2lT aN< ^ 2 ^ b pxG{x, 1/x 2 ) 



(49) 



which is exactly formula (21b). Thus one indeed is able to calculate the produced gluon spectrum from the 
"classical" field given in (37) and (38) and this accounts for the final state interactions which are present in 
covariant gauge. 

3 The produced gluon distribution in nucleon-nucleus collisions 

In this section we shall calculate the spectrum of gluons produced in the scattering of a nucleon off a large 
nucleus. If the nucleus is right-moving and the nucleon left-moving then we suppose the centrally produced 
gluon momentum, £, is such that £- is much less than the minus component of the nucleon momentum and 
that £ + is sufficiently less than the momentum of a nucleon in the nucleus so that the gluon is coherent over 
the longitudinal extent of the nucleus. We begin our discussion with gluon production in nucleon-nucleon 
collisions. 

3.1 Soft gluon production in nucleon-nucleon collisions 

We choose A- = gauge and a center of mass system to describe gluon production. If we choose the 
outgoing gluon polarization such that rj ■ e = e + = then we need only consider the graph shown in Fig. 8 
for the gluon production amplitude. The terms occurring at the 3-gluon vertex in that graph are 

V W = k _ {£ \ (£ - k)^ef[-g^{2k - £)p + gp a {k + l) 7 - g^(2£ - k) a ] (50) 
where pi+ and p 2 - are the large components of pi and P2 respectively. It is straightforward to check that 

v {x) = \lh {x) ■(£-£)- if - fc 2 k (A) -l] (si) 



and that 



Y V w V M* = *k 2 (£-V)\ 



There is an fj p at the p— vertex in Fig. 8 while there is a k^r at the a— vertex. Writing 

= k a - r\ a k^ - fj a k + (53) 

and using the fact that k a gives zero by current conservation while rf a gives a very small result one can 
replace k„ by —rj a k^. The fc_ will cancel the denominator in (51) (and in (52)) so that in each case the 
vertices p and a are multiplied by the appropriate factors for the cross section to be given by the product of 
the unintegrated gluon distributions of nucleons p\ and pi- The result is]l4]|-[|T9[. 

da 4aN c f 2 dxtGjxutf) dx 2 G(x 2l (£ - k) 2 ) 



<P£dy (N 2 -l)fJ dk 2 d{£-kf 

where xix 2 s — f 2 and y is the rapidity of the produced gluon. Of course (54) is difficult to take literally. If 
f 2 is not large there is no hard scale and the unintegrated gluon distributions do not have much meaning. 
On the other hand if £ 2 is large then the cross section is better expressed as the product of two gluon 
distributions times a hard scattering cross section. Nevertheless, (54) is useful as a normalization and for 
comparison with the equations we are now going to derive for gluon production in nucleon-nucleus collisions. 
It also follows in a BFKL approximation pO, UB. 



11 



3.2 Central region gluon production in nucleon-nucleus collisions 

It is not too difficult to calculate the cross section for producing a gluon of transverse momentum I in the 
collision of a nucleon with a large nucleus. It is convenient to view the calculation in the rest system of the 
nucleus. We suppose the longitudinal momentum of the gluon is large enough so that the gluon is coherent 
over the length of the nucleus Jl3[. We begin by considering gluon radiation when a quark scatters on a 
nucleus. It will then be a simple matter to extend the result to an incoming nucleon. It is convenient to 
work in fj ■ A = A- = gauge and with fj • e = where we choose the incoming quark to have a large minus 
component of momentum. Suppose the quark reaches the nucleon at light-cone time r = = 0. Then we 
separate the gluon radiation cross section into four terms depending on the time T\ that the gluon is emitted 
from the quark in the amplitude and on the time r 2 that the gluon is emitted from the quark in the complex 
conjugate amplitude. The four cases are: 



(a) ti > 0, r 2 > (55a) 

(b) n < 0, t 2 > (556) 

(c) ti > 0, t 2 < (55c) 

(d) ti < 0, r 2 < (55d) 

Since the coherence time of the gluon is assumed to be much bigger than the radius of the nucleus we can 
assume that the quark, or the quark-gluon system, passes over the nucleus instantaneously compared to 
the magnitudes of the times of emission t\ and r 2 . Also once we separate contributions into definite times 
we are in effect dealing with light-cone perturbation theory rather than with Feynman graphs though in 
many graphs the orderings (55) simply tell whether the gluon is emitted before or after the quark, or quark- 
gluon system, interacts with the nucleus. The cross sections coming from the regions (55) are given by the 
expressions 

^1 = I f d 2 bd 2 xid 2 xo — aCF - 1 IgM fe-g.) (56) 
dHdy 7r J aMXiax2 ^ 7T x\x 2 [ ' 



d<? (b+C) = - I [tf^ 1 aC F £l-£2 

d 2 £dy ir J 1 2 An 2 n x\x\ 



, x\v\J B? - b 2 , , xlv^R 2 - b 2 A Mt _ x s .„_. 
expt- - 1 2A =-] + cxp[- ~ 2 V 2A =-] e^fe a) (57) 



^zfez^EZ,^ (58 , 

In the above we take the transverse coordinate of the quark to be o while x^ and x 2 are the transverse 
coordinates of the gluon in the amplitude and complex conjugate amplitude respectively. We use a shorthand 
notation where x 2 v means x 2 v(x 2 ) with v, A, R and b as previously used in Sec. 2. Let us now see how (56) 
to (58) come about. 

In (56) the factor ^f- ~^i is just the product of the coordinate space gluon emission amplitude 

times the complex conjugate amplitude. What may seem surprising in (56) is that there is no trace of the 
nucleus! But this is easy to understand. What passes through the nucleus is a high energy quark. The quark 
does not change its transverse coordinate nor does it lose a significant amount of energy as it passes through 
the nucleus. (Technically, elastic quark-nucleon scatterings and inelastic quark-nucleon scatterings, where 
the nucleon breaks up, cancel.) What is unusual in (56) is that we are attributing a contribution to the 
cross section to the disconnected graph where the quark freely passes through the nucleus. The Feynman 
diagram contribution to this process would be zero by energy conservation. However, (56) corresponds to 
T\ and Ti greater than zero so we do not yet have energy conservation. There are other terms in (57) and 
(58) which can be viewed as cancelling this "gluon radiation without interaction" term. Indeed these are the 
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^-independent terms of (57) and (58). However, the grouping of terms as given in (56)-(58) is convenient for 
our purposes. 

Now turn to (58). The exponential term in (61) reflects multiple scattering of the gluon as it passes 
through a length Az — 2\/ R 2 — b 2 of nuclear matter. This factor is identical to that given by (20) with 
z - z taken to be 2^R 2 -b 2 . In this case it is the quark-gluon system which passes through the nucleus 
both in the amplitude and in the complex conjugate amplitude. Interactions of the quark with the nucleons 
in the nucleus cancel between real and virtual (production and elastic scattering) terms. 

The expression (57) corresponds to a quark-gluon system passing through the nucleus in the amplitude 
and a quark passing through the nucleus in the complex conjugate amplitude along with a term where 
the amplitude and complex conjugate amplitude terms are exchanged. In Appendix A we outline how the 
exponential factors, involving x 2 and x 2 , come about. 

It is straightforward to evaluate (56)- (58) in the approximation of neglecting the x-dependence in v. One 
finds 

da (a ) f aC F 1 

d 2 b (59) 



d 2 ldy J ir 2 I 2 
da^ f aC F [1 - 



d 2 £dy J it 2 £ 2 



l -d 2 b (60) 



= / d z b 



d 2 Uy J it 2 <£ 2 > 



< £ 2 > L 2 £ 2 -I 2 



whereflll 



< I 2 >=< £ 2 {b) >= 2 J IP - b 2 H< [ >] . (62) 

A 

L is an infrared cutoff and V(n, z) is the incomplete T- function. In arriving at (61) we have used 



o 



'My)7j^ = 2 ^(*-2>)| ^ 



and 

r-2n 



1 2tt 

\x + y| 2 \x 2 — y 2 

with x and y in (63) and (64) related to x x and x 2 in (58) by The second plus third terms 

in the brackets on the right-hand side of (61) are non-logarithmic when i 2 <C< £ > 2 . When I 2 ^>< £ 2 > this 
term exactly cancels similar terms in (59) and (60). Thus in the logarithmic approximation we keep only 
the first term on the right-hand side of (61) and we understand that (59) holds only for £ 2 £ 2 > . 
We can now go to the case of an incident nucleon by making the replacements 

^^^G(x,£ 2 ). (65a) 

( ^£n[<£ 2 >L 2 \^xG{x,£ 2 ). (656) 
In the logarithmic approximation the (b) + (c) contribution can be neglected and one finds 



da = i / d 2 b 



d 2 £dy tt 



f) exp(-^|U) 
— xG(x, £ 2 ) + xG(x, < £ 2 » J 



(66) 
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with the understanding that the first term on the right-hand side of (66) can only be used when £ 2 <C< £ 2 > ■ 
We note that 

T v =JJk' ( ' t = 2 J d ' bxG ^ <e>) (67) 

There is a simple physical interpretation of (66) and (67). When an incident nucleon passes through a 
nucleus all the gluons in the nucleon's wavefunction having £ 2 < < £ 2 > are freed during the collision and 
the final state interactions broaden their distribution. This is the second term on the right-hand side of 
(66) and one half of (67). The nucleon remnants emerging from the nucleus have lost their gluon cloud. In 
rebuilding that cloud there is further gluon emission, but this gluon emission has no transverse momentum 
broadening. This is the first term on the right-hand side of (66) and the other half of (67). 



3.3 Gluon production in nucleon-nucleus collisions without final state interac- 
tions 

In this section we describe exactly the same process as in the previous section. However, now we choose 
A + = gauge with e_ = 0. Our object is to show that, as for our deep inelastic scattering process of Sec. 2, 
the gluon production cross section can be described completely in terms of the interaction of the "classical" 
fields of the incoming nucleon and of the incoming nucleus without any final state interactions. 

It is relatively straightforward to derive a general formula for d i[ d in terms of classical fields associated 
with the colliding nucleon and nucleus. To see how this comes about we discuss in detail the scattering of 
a left-moving nucleon on two right-moving nucleons. We begin with the graphs shown in Fig. 9 where we 
suppose that nucleons 1 and 2 are right-moving and that the x_— coordinate of nucleon 1 is less than that of 
nucleon 2, that is, nucleon 1 is ahead of nucleon 2. We assume these nucleons are separated by a longitudinal 
distance large compared to 1/A in their common rest system. Then the interaction in Fig. 9a must be a gauge 
term with the effective part of the k-propagator being p^rr^] ■ The phase factor e lk +(. x 2~ x i)- dictates that 
one distort the contour in the upper-half plane with the propagator becoming ^gfS(k + )k a . Similarly for the 
graph in Fig. 9b one can again distort the propagator in the upper half fc+— plane getting an identical result 
for the fc— propagator. One now can use the STW identities to arrive at a contribution which is illustrated 
in Fig. 10 where the dotted line brings color to the (£ — k — r) line as well as transverse momentum, but it 
brings no longitudinal momentum and no momentum dependence at its vertex. This can be interpreted 
as a color rotation and transverse momentum convolution of the field of nucleon 1 by the gauge field from 
nucleon 2 followed by a scattering off the left-moving nucleon. 

Now consider the graphs shown in Fig.ll. We again distort the fc+ — contour of integration into the 
upper half plane. In Fig. 11a there is also the light-cone denominators from the (fc + r)— line which must 
be considered. The term having a singularity in the upper-half plane is proportional to jj^j-^j^- It is 
straightforward to check, using the STW identities, that this term gives zero because of our choice fj ■ e = 0. 
Thus in graph(a) there is only the k ^" ie singularity of the fc— line which is evaluated by the STW identities. 

The graph shown in Fig. lib has two contributions, one being the k ^ ie term from the fc— propagator and 
the second being a singularity in the nucleon-gluon scattering process (P) + (— r — k) — > (P — r) + (— fc). 
This second term corresponds to a single gluon exchange final state interaction of the left- moving nucleon 
remnants with nucleon 2 and it cancels with corresponding two-gluon exchange terms in the amplitude and 
complex conjugate amplitude. The result of the terms in Fig.ll is illustrated in Fig. 12 and can be interpreted 
as a color rotation, and convolution, of the field from the left-moving nucleon by nucleon 2 followed by a 
scattering off nucleon 1. 

Now we consider the graphs shown in Fig. 13. We shall evaluate these graphs by distorting the fc + — contour 
into the lower half plane. Since we now distort in a direction opposite to that natural for the factor 
e ik + (x 2 -x 1 )- jj. j g inypoj-tant to check that there is sufficient convergence in the fc+ — plane so that the 
exponential factor is not relevant. For example, in graph(a) one has fc + — denominators [fc+ — ie] _1 , [(fc + r — 

t)+ + ie]^ 1 and [(fc — £) + + ie]^ 1 multiplying fc^ and (fc + r — £)^ as well as a possible light-cone 

denominator from the {£— fc)— line. (Note, however, that the ^'-e)'++ie term m the (^ — fc)— propagator gives a 
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small contribution because of our choice of polarization,^- e = 0.) Thus there is sufficiently strong convergence 
for the k + — integral unless (possibly) there is a factor of fc+ coming from each of the 3-gluon vertices. At the 
r^«7 vertex the only possibility leads to a factor of ij p . (Factors of fj a or fjj give zero because they multiply 

or e 7 respectively.) At the T\ pu vertex a factor of fj v leads to an extra denominator \(k — as does a 

factor of fj\ since (k + r — £)pg pv f} v = 0. Thus there is sufficient convergence to distort the k + — contour in the 
lower half plane where one picks up only the pole \{k + 1 — £) + + ie]^ 1 . After picking up this pole one writes 
(k + r — £)p = (k + r — £) p — r]p(k + r — £)_ — fip(k + r — l) + and using (k + r — £) + = and (k + r — £)_ w one 
again can use the STW identities as was done for the graphs of Fig.ll. Here one obtains a result illustrated 
in Fig. 14 along with a rescattering of the remnant of the left-moving nucleon which will cancel as discussed 
above. Hence, we may view the graphs of Fig. 13 as a color rotation, and transverse momentum convolution 
of the field from the left-moving nucleon by nucleon 1 followed by a scattering off nucleon 2. 

Finally, there is the graph shown in Fig. 15. Here we distort the contour in the upper half plane to pick 
up the pole at k + = with the factor k a acting at the 4-gluon vertex after the k+ — integration has been 
done. 

In order to see that the graphs of Figs. 9, 11, 13 and 15 give the same result as our previous calculation 
summarized in (58) it is convenient to evaluate the graphs of Fig. 13 in a different way. If we evaluate these 
graphs by distorting the k + — contour into the upper half fc + — plane poles at k + = and at {£ — k) + = 

are encountered. (There is also the term rj^jF^zrii comm g fr° m the {£ — k)— line propagator which gives 
an additional singularity, however this term leads to a small result since the (k — £) u term either cancels 
between the a and b parts of Fig. 13 by STW or the (k — £) v eliminates the (£ — k — r)— propagator leaving 
no singularities whatever in the lower half k + plane.) When the singularity at k + = is taken one has a 
resulting k a which when combined with the results given by the graphs in Figs. 10, 12 and 15 gives zero from 
the STW identities. The singularity at (£ — k) + — ^p gives on-shell propagation of the {£ — A:)— line with 
a rescattering off nucleon 2 after production by nucleon 1. Thus we recover the result of our calculation in 
Sec. 3. 2 when that calculation is restricted to the interaction of just two nuclcons in the nucleus. 

However, one can also describe the scattering in terms of the "classical" A+ = light-cone gauge fields 
of the left-moving nucleon and of the 2 right-moving nuclcons. In this case the left-moving nucleon has a 
classical field color rotated by nucleons 1 and 2 as indicated in Figs. 12 and 14. The right-moving nucleons 
have their Weizsacker- Williams fields where the field of nucleon 1 is rotated by that of nucleon 2. In every 
case the nucleons with larger values of x_ rotate the fields coming from smaller values of x_ according to 
(37) and (38). The gluon production amplitude is then given by the three and four gluon interactions of 
these incoming "classical" fields as described by the graphs in Figs. 10, 12, 14 and 15. In Appendix B we 
shall outline the technical argument for this result when an arbitrary number of nucleons in the nucleus are 
involved in the interaction. 

We now put into formulas what we described above for gluon production in a nucleon- nucleus collision 
calculated in A + = gauge, and where the nucleon is left-moving and the nucleus is right-moving. We 
denote by A^(k,k + ) the field of the nucleus given by (37-40) and (42). The "classical" field of the nucleon 
is denoted by 

A'^k,k-,b) = A'(k,k-,b)r lli (68) 

with 

A'(x,x+,b) = -S{x)T a S- 1 {x)5{x + )ln[\x-b\n]p a N (69) 
where S(x) is the same as in (38) but with X- = — oo, and where 

< p a N p a N >= j^jQ 2 ^ xG ^ Q 2 ) ( 7 °) 

indicates an average over internal nucleon structure. We have used current conservation at the nucleon 
source to eliminate the light-cone gauge denominator in A' p much as was done in the discussion following 
(53) and which led to (54), but now with the roles of + components and — components exchanged. In order 
to write a compact formula for d i% d it is convenient to define a total gluon field A 1 " 1 (x) as 
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A™(x) = A^x) + A'^x,b) + A^{x) (71) 

where A 1 - and A' are the fields given in (37) and (68), respectively, while A^~ ee {x) is a free quantized gluon 
field normalized according to 

(A) / n\ i£-x 

^ A 'r^ = 7pfc r ° (72) 

with \£Xa) being a single gluon state having momentum £, polarization A and color a. A tot depends on the 
impact parameter of the collision through A'. Then the gluon production cross section is 

, „ 2 N 2 C -1 

° =LUe d%<J2 E (0|S|*Aa)(*Aa|S| 0) > (73) 

y J A=l a=l 



with 



S =J- I d 4 x TrF^ (x)F flv (x) (74) 



2 

and where 

F\iv = dpA? - 8 V A^ - ig[A^\ Ai ot ] . (75) 

The <> in (73) indicate an average over nuclear structure as indicated in (39) as well as an average over 
nucleon structure as given by (70). In evaluating (73) one assumes that the field A'J ee is much smaller than 
A'^ and Ah. Therefore A^ ree is included only once in the 3 and 4 gluon vertices in S. At the same time each 
of the other two fields has to be included at least once to allow the production of the gluon. Eq.73 is a 
compact, and elegant, formula for gluon production in nucleon-nucleus collisions. 



16 



Appendix A 



In this appendix we shall show how the exponential factors in (57) come about. We consider the case where 
Ti < and T2 > 0. We consider the last, the latest in time, interaction with a nucleon in the nucleus. 
In the amplitude the interaction is with a quark-gluon system with the quark at transverse coordinate 
and the gluon at x 1 . In the complex conjugate amplitude the interaction is only with the quark, again at 
transverse coordinate 0. The possible interaction terms, both real and virtual, are shown in Fig. 16. The 
factors associated with each graph are listed below as 



(«) = 


-C F /2 V(Q)/N c 


(b) = 


N c /2 V{x x )/N c 


(c) = - 


-1/2N C V(Q)/N C 


(d) = 


-N c /2 V(Q)/N C 


(e) = 


N c /2 V{x 1 )/N c 


(/) = 


-C F /2 V{Q)/N c 



where V is given in (17). Adding all the terms one finds 

(a) + (6) ■ ■ ■ +(/) = -\V(P) = -\x\vix,). (A.2) 

Multiplying the right-hand side of (A.2) by per = 1/A to take account of the density of nuclcons and the 
cross section we obtain the same type of factor as appears on the right-hand side of (14). Our convention is 
to use a, and A, for gluon scattering on nucleons. The factor of 1/N C on the right-hand side of the terms in 
(A.l) converts the result from gluons to the appropriate partons involved in the interaction. 



Appendix B 

In this appendix we outline an argument for the absence of final state interactions in gluon production 
in nucleon-nucleus collisions. The major barrier in generalizing the discussion given in Sec. 3. 3, for the 
absence of final state interactions when a nucleon scatters on two nuclcons, is the presence of light-cone 
gauge denominators of the type ^ k _^ +ie , along the gluon line carrying the large minus-component of the 
momentum, which hinder the distortion of k + — contours into the lower half plane. In Sec. 3. 3, for the graphs 
in Fig. 13 this term was eliminated by our choice of polarization vector, fj ■ e = for the produced gluon. In 
the general case the troublesome graphs are of the type shown in Fig. 17. The graph in Fig. 17a is such that 
cutting the (£— fc)-line separates the overall graph into two parts while the graph in Fig. 17b is an example of 
a graph where cutting the (£ — fc)-line does not separate the graph into two parts. There is not any particular 
ordering of the x_ — positions of the nucleons belonging to A and B. That is, some nucleons in A may 
have larger values than some of the nucleons in B. We begin with the graphs of the type shown in 
Fig. 17a. For our purposes the dangerous term is the (2-i)++ie P art °^ tnc — ^) — propagator. However, 
once all hookings of the (t — k) — line in B are included, at a given order of the coupling and with a given 
set of nucleons, the net result must be zero from the STW identities. 

Now turn to graphs of the type shown in Fig. 17b where one or more lines, in addition to the (£— k)— line 
connect A and B. A particular example of such a graph is shown in Fig. 18, where nucleon 1 and nucleon 
3 have inelastic reactions while nucleon 2 has an elastic reaction. The potentially dangerous problem with 
graphs of this type is that the application of the STW identities does not directly give zero since the k\ line 
carries color. Indeed, applying the STW identities for the term in the (£—k)— propagator leads to 

a term like that indicated in Fig. 19 where the dashed line carries the propagator factors ^ k __^yi +it ^ k Zf)l +it 
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and inserts a momentum l~ k on the k\ — line as well as rotating the color factors of the k\— line. However, 
this term is small because the scattering with nucleon 2 is the same as the high energy scattering gluon 
(£ — k — fci)+ nucleon — > gluon (I — ^2)+ nucleon and the high energy limit of such a scattering is zero in 
our approximation of no evolution within nucleon 2. 

Thus, the gluon line carrying the large minus component of the momentum £, the hard gluon line, will 
never give singularities as one distorts k + ~ momenta into the lower half plane. This is all that is needed 
to arrive at our main result, described in Sec. 3. 3 and summarized in (73)-(75). To see this, begin with the 
nucleon having the largest values of X-. A gluon, g, from this nucleon either attaches directly to the hard 
gluon line or to a gluon coming from another nucleon or, perhaps, directly to another nucleon. If the gluon, 
<?, does not attach directly to the hard gluon one is in the case considered in Rcf.ll where it was shown 
that g rotates the fields coming from nucleons with smaller values of x_. If the gluon, g, attaches to the 
left-moving nucleon, or to the hard gluon before the hard gluon has any previous interactions, one distorts 
the k + of g into the upper half plane and uses the STW identities. If g attaches to the hard gluon after the 
hard gluon has interaction with other nucleons one routes the momentum of g through the hard gluon and 
back into a nucleon at the first available interaction, as illustrated for example in Fig. 18. In this case one 
distorts the — contour in the lower half plane and uses the STW identities not for g but for the gluons 
carrying the momentum k returning from the hard gluon. This set of operators leads to the fields and 
A'^ in (71) and a cross section given by (73). 
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